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Abstract 

We use two-dimensional Poissonian random lattices of Voronoi/ 
Delaunay type to study the effect of quenched coordination number 
randomness on the nature of the phase transition in the eight-state 
Potts model, which is of first order on regular lattices. ^From extensive 
Monte Carlo simulations we obtain strong evidence that the phase 
transition remains first order for this type of quenched randomness. 
Our result is in striking contrast to a recent Monte Carlo study of 
quenched bond randomness for which the order of the phase transition 
changes from first to second order. 
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1 Introduction 


For systems exhibiting a continuous phase transition in the pure case it is 
well known that the influence of quenched random disorder can modify the 
critical behaviour, leading to new universality classes or even eliminating 
the phase transition altogether [1], Also for systems undergoing a first- 
order phase transition the effect can be very dramatic. Phenomenological 
renormalization-group arguments suggest that the addition of quenched ran¬ 
domness can smoothen the transition completely and induce instead a con¬ 
tinuous phase transition [2], For a certain type of quenched bond-disorder 
in the two-dimensional g-state Potts model with q = 8 this prediction has 
recently been confirmed by extensive Monte Carlo simulations [3]. While the 
pure model is exactly known to exhibit a quite strong first-order phase tran¬ 
sition [4], the simulations with quenched bond-disorder gave clear evidence 
for a continuous phase transition. From careful finite-size scaling analyses 
the phase transition was identified to belong to the Ising model universality 
class. 

In this note we report Monte Carlo simulations of the same model sub¬ 
ject to a different kind of quenched disorder. Instead of using a square lattice 
with uniform coordination number (=4) and drawing the coupling strengths 
randomly from two different values as in Ref. [3], we consider Poissonian ran¬ 
dom lattices where the coordination numbers vary locally between 3 and oo 
and the coupling strengths are uniform. The random lattices are constructed 
according to the Voronoi/Delaunay prescription [5] for toroidal topology, i.e., 
with periodic boundary conditions. 


2 Model and simulation 

We used the standard definition of the g-state Potts model, 

Zpotts = = ( 1 ) 

Wi} (*J> 


where (3 = JjUsT is the inverse temperature in natural units, and {ij) de¬ 
notes the nearest-neighbor bonds of random lattices with V = 250, 500, 750, 
1000, 2000, and 3000 sites. For each lattice size we generated 20 indepen¬ 
dent replica and performed long simulations of the 8-state model near the 
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transition point at /3 = 0.826, 0.830, 0.830, 0.830, 0.832, and 0.833, respec¬ 
tively, using the single-cluster update algorithm [6]. After thermalization we 
recorded 1 000 000 measurements (taken after 1, 1, 1, 1, 2, 4 clusters had been 
flipped) of the energy and the magnetization M = (gmax{ni} — 14)/(g —1) 
in a time-series file, where rii < V denotes the number of spins of “orienta¬ 
tion” i = 1,. . . in one lattice configuration. Obviously it is sufficient to 
store the integers Vjq < max{ni} < V. From this data it is straightforward 
to compute all quantities of interest as a function of temperature by stan¬ 
dard reweighting procedures [7]. The corresponding quantities per site are 
denoted in the following hy e = EjV and m = MjV. 

More precisely we used reweighting to compute, e.g., the specific heat 
= /3^V ((e^) — (e)^) for each replica labeled by the superindex (f), and 
then performed the replica average C'(/3) = [C'^*^(/3)] = (1/20)'^^° 
denoted by the square brackets. To perform the replica average at the level 
of the (7^*^ (and not at the level of energy moments) is motivated by the 
general rule that quenched averages should be performed at the level of the 
free energy and not the partition function [8]. Finally, we determined the 
maximum, Cmax = for each lattice size and studied the finite- 

size scaling (FSS) behaviour of Cmax and The error bars on the two 

quantities entering the FSS analysis are estimated by jack-kniving [9] over 
the 20 replicas. This takes into account the statistical errors on the estimates 
of each as well as the fluctuations among the different caused 

by the quenched randomness. The analysis of the magnetic susceptibility, 
x(/3) = 13V {[{vn?) — {rnY\) proceeds exactly along the same lines, yielding 
Xmax and . 

In the case of the (energetic) Binder parameter, usually defined on reg¬ 
ular lattices as 5(/3) = 1 — (e^)/3(e^)^, the proper definition of the replica 
average is less clear to us. In order to study this problem we have therefore 
computed the following three definitions which differ only by the replica aver¬ 
aging procedure: 5i(/3) = 1 - [(e^)/3(e2)2], 52(/3) = 1 - [(e^)] /3 [(e^)^], and 
= 1 — [(e^)] /3 [(e^)]^. While in spin glass simulations [10] usually the 
analogue of B 3 (with e replaced by the overlap) is used, for a random bond 
Ising chain [11] a better scaling behaviour was observed for the analogue of 
Bi (with e replaced by m). 
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3 Results 


Already a first qualitative inspection of our data gave a clear indication that 
the first-order nature of the phase transition on regular lattices (square, tri¬ 
angular, . . .) persists on quenched random lattices. To make this statement 
more quantitative let us first consider the FSS of the specific-heat and suscep¬ 
tibility maxima. If the hypothesis of a first-order phase transition is correct, 
we expect for large system sizes an asymptotic FSS behaviour of the form 
[12-14] 

C*max = O-C -\- bcV + . . . , ( 2 ) 

and 

%max bx^ T • • • • (3) 

Our data shown in Fig. 1 are clearly consistent with the Ansatz (2), (3). The 
least-square fits yield ac = 23.3(2.0), be = 0.0659(30), with a goodness-of-fit 
parameter Q = 0.16 (corresponding to a chi-square per degree of freedom of 
1.7), and = -0.70(43), 6 ^ = 0.0629(13), with Q = 0.45. 

In Fig. 2 we show the scaling of the Binder parameter minima which is 
expected to be of the form 

-Si,min = 0-Bi + bBijV + . . . . (4) 

Since the data for Bi and B 2 are almost indistinguishable, we have only shown 
Bi. Again the data confirms the hypothesis of a first-order phase transition, 
and from the fits we obtain = 0.6240(20), bs^ = —18.8(1.4), Q = 0.17, 
us, = 0.6236(22), be, = -18.5(1.4), Q = 0.47, and ae, = 0.61125(68), 
6^3 = —16.45(71), Q = 0.55. Notice the much higher accuracy of B 3 . 

The locations of the extrema of C{j3), x(/3)) ^'iid Bi[/3) define pseudo¬ 
transition points which, at a first-order phase transition, should scale ac¬ 
cording to 

/^Cmax = /3o + CejV + . . . , (5) 

etc., where /3o is the infinite volume transition point. Our data and the 
corresponding fits through all data points are shown in Fig. 3. The resulting 
estimates for /3o are 0.83360(14) from Cmax {Q = 0.51), 0.83365(14) from 
Xmax {Q = 0.47), and 0.83362(13) from 53 ,mm {Q = 0.23). On the scale 
of Fig. 3 the data points for 5i,mm and 52 ,mm could hardly be disentangled 
from 53 ,mm and are therefore omitted. The results for /3o are 0.83371(14) 
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from {Q = 0.40) and 0.83350(13) from 52,mm {Q = 0.25). Taking the 

average of the different estimates we finally obtain 

/3o = 0.83362 ± 0.00013. (6) 

Notice that this value is very close to the exactly known transition point of 
the 8-state Potts model on a triangular lattice = 0.85666 . . .) [4]. 

Finally we had a closer look at the replica fluctuations. In Fig. 4(a) 
we show the curves for all 20 replica as well as the resulting replica 

average C'(/3) = [C'0)(/3)]. We see that all curves look very similar but are 
displaced by a constant amount in j3. This is illustrated in Fig. 4(b) where 
we plot the same data vs /3 —where (3 ^1^ is the location of the specific- 
heat maximum for the i’th replica. Here all curves fall almost on top of 
each other and define quite nicely a kind of master curve. This suggests 
that the main effect of the randomness in the coordination numbers is to 
define a temperature reference point that is fluctuating among the different 
replica. This makes it plausible that coordination number randomness is not 
sufficiently strong to change the nature of the phase transition. 


4 Conclusions 

Summarizing, we have obtained clear numerical evidence for a first-order 
phase transition in the 8-state Potts model on quenched random lattices of 
Voronoi/Delaunay type. We can savely exclude the possibility of a cross-over 
to a continuous transition as was observed for a certain type of quenched bond 
randomness on square lattices [3]. 

In this brief note we have confined ourselves to FSS analyses of standard 
observables. It would be interesting to extend the analysis to quantities that 
are directly related to the probability distributions of the energy or mag¬ 
netization, such as the interface tension and the “ratio-of-weight” definition 
of pseudo-transition points (which gave very precise estimates for regular 
square lattices) [14]. A quite elaborate study in this direction based on a 
much larger set of 128 replica will be published elsewhere [15]. 


4 



Acknowledgements 

W. J. would like to thank D.P. Landau for useful discussions and communicat¬ 
ing the results of Ref.[3] prior to publication. W.J. and R.V. were supported 
in part by NATO grant CRG940135. W.J. thanks the Deutsche Forschungs- 
gemeinschaft for a Heisenberg fellowship, and also acknowledges support in 
part by EC grant ERBCHRXCT930343. The Monte Carlo simulations were 
performed on clusters of fast RISC workstations in Barcelona and Mainz. 


References 

[1] A.B. Harris, J. Phys. C7 (1974) 1671; Y. Imry and S.-k. Ma, Phys. 
Rev. Lett. 35 (1975) 1388; A. Aharony, Phys. Rev. B18 (1978) 3318; 
C. Crinstein and S.-k. Ma, Phys. Rev. Lett. 49 (1982) 684; A.N. Berker, 
Phys. Rev. B29 (1984) 5293. 

[2] A.N. Berker, Phys. Rev. B29 (1984) 5293; K. Hui and A.N. Berker, 
Phys. Rev. Lett. 62 (1989) 2507; ihid. 63 (1989) 2433(E); A. Aizenman 
and J. Wehr, Phys. Rev. Lett. 62 (1989) 2503; A.N. Berker and K. Hui, 
in Science and Technology of Nanostructured Magnetic Materials, eds. 
C.C. Hadjipanayis, C. Prinz, and L. Paretti (Plenum, New York, 1991). 

[3] S. Chen, A.M. Eerrenberg, and D.P. Landau, Phys. Rev. Lett. 69 (1992) 
1213 (1992); and UGA preprint (1994). 

[4] E.Y. Wu, Rev. Mod. Phys. 54 (1982) 235; 55 (1983) 315(E). 

[5] J.L. Meijering, Philips Res. Rep. 8 (1953) 270; R. Collins, Proc. Phys. 
Soc. 1 (1968) 1461; and in Phase Transitions and Critical Phenom¬ 
ena, Vol. 2, eds. C. Domb and M.S. Green (Academic Press, London, 
1972), p. 271; J.M. Ziman, Models of Disorder (Cambridge University 
Press, Cambridge, 1976); C. Itzykson and J.-M. Drouffe, Statistical Field 
Theory (Cambridge University Press, Cambridge, 1989), Vol. 2; N.H. 
Christ, R. Eriedberg, T.D. Lee, Nucl. Phys. B202 (1982) 89; Nucl. Phys. 
B210 [PS6] (1982) 310, 337; M. Tanemura, T. Ogawa, and N. Ogita, 
J. Comp. Phys. 51 (1983) 191; C. Itzykson, Fields on Random Lattices, 
in Progress in Gauge Field Theories, proceed. Cargese Summer School, 


5 



ed. G.’t Hooft (Plenum Press, New York, 1983); R. Friedberg and H.-C. 
Ren, Nucl. Phys. B235 [FSll] (1984) 310; H.-C. Ren, Nucl. Phys. B235 
[FSll] (1984) 321. For a recent study of the Ising model on 2D random 
lattices of this type see, W. Janke, M. Katoot, and R. Villanova, Phys. 
Rev. B49 (1994) 9644. 

[6] R.H. Swendsen and J.S. Wang, Phys. Rev. Lett. 58 (1987) 86; U. Wolff, 
Phys. Rev. Lett. 62 (1989) 361. 

[7] A.M. Ferrenberg and R.H. Swendsen, Phys. Rev. Lett. 61 (1988) 2635; 
ihid. 63 (1989) 1658(E). 

[8] K. Binder and A.P. Young, Rev. Mod. Phys. 58 (1986) 801. 

[9] R.G. Miller, Biometrika61 (1974) 1; B. Efron, The Jackknife, the Boot¬ 
strap and other Resampling Plans (SIAM, Philadelphia, PA, 1982). 

[10] E. Marinari, G. Paris!, and E. Ritort, Roma preprint (1993), cond- 
mat/9310041. 

[11] A. Crisanti and H. Rieger, J. Stat. Phys. 77 (1994) 1087. 

[12] C. Borgs and R. Kotecky, J. Stat. Phys. 61 (1990) 79; Phys. Rev. Lett. 
68 (1992) 1734; C. Borgs, R. Kotecky, and S. Miracle-Sole, J. Stat. Phys. 
62 (1991) 529; J. Lee and J.M. Kosterlitz, Phys. Rev. B30 (1991) 3265; 
W. Janke, Phys. Rev. B47 (1993) 14757. 

[13] Eor a general overview see, e.g., K. Binder, Rep. Prog. Phys. 50 (1987) 
783; or the articles in Dynamics of First Order Phase Transitions, eds. 
H.J. Herrmann, W. Janke, and E. Karsch (World Scientific, Singapore, 
1992). 

[14] See, e.g., W. Janke, in Computer Simulations in Condensed Matter 
Physics VII, eds. D.P. Landau, K.K. Mon, and H.B. Schiittler (Springer 
Verlag, Heidelberg, Berlin, 1994), p. 29; and references therein. 

[15] W. Janke and R. Villanova, in preparation. 


6 









Figure 4: (a) The specific heat for each of the 20 replica as a function 
of inverse temperature, and the replica average C = [(7^*^] (marked by filled 

(b) The same data plotted vs j3 — where denotes the maximum 

location for the i ’th replica. This shows that the main effect of the randomness 
in the coordination numbers can be parametrized by a random temperature 



offset. 
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